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ABSTRACT 

We present Chandra X-ray Observatory archival observations of the supernova rem- 
nant 1E0102. 2-7219, a young Oxygen-rich remnant in the Small Magellanic Cloud. 
Combining 28 ObsIDs for 324 ks of total exposure time, we present an ACIS image 
with an unprecedented signal-to-noise ratio (mean S/N ~ \/S ~ 6; maximum S/N > 
35) . We search within the remnant, using the source detection software wavdetect, 
for point sources which may indicate a compact object. Despite finding numerous de- 
tections of high significance in both broad and narrow band images of the remnant, we 
are unable to satisfactorily distinguish whether these detections correspond to emission 
from a compact object. We also present upper limits to the luminosity of an obscured 
compact stellar object which were derived from an analysis of spectra extracted from 
the high signal-to-noise image. We are able to further constrain the characteristics of a 
potential neutron star for this remnant with the results of the analysis presented here, 
though we cannot confirm the existence of such an object for this remnant. 



1. Introduction 

The sup ernova remnant IE 0102 .2-7219 (SNR 1E0102. 2-7219, hereafter, "E0102") was first 
identified by ISeward Sz Mitchelll ( 198ll ) in a survey of the Small Magellanic Cloud (SMC) with the 
Einstein X-Ray Obs ervatory. While this oxygen-rich remnant resides at a distance of 59 kpc 
(jvan den Berghlll999l ). it is far less extincted than many Galactic remnants due to its high galactic 
latitude (/ ~ -45°). Kinematic studies in the X-ray and optical have determined E0102 to be 



a rel a .tively young remnant , with an expansion age between 1000 and 2000 years (JHughes et al. 



2000 



Finkelstein et al.l l2006l . respectively). For these reasons , E0 102 has been extensively studied 
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E0102 was likely generated by a Type lb or Ic core-collapse supernova event (JBlair et al.ll200d ). 
As such, the remnant E0102 would be expected to host sor ne variety of compa ct o bject (i.e., stellar 
mass black hole or degenerate stellar object). Furthermore. lAmv &: Balll (| 19931 ) and lFinkelstein et al. 
(|2006i ) have considered the possibility that such a compact object may be influencing the complex 
morphology of this remnant. Though this remnant has been well studied, few groups ha ve provided 
observ ational constraints to the presence of a compact object for E0102. In fact, only iGaetz et al. 
([2000|) have presented an upper flux limit on a Crab-like pulsar for this remnant. In this research, 
we seek to constrain the X-ray luminosity of this unidentifled compact object. 

E0102 is a calibration target for the Chandra X-Ray Observatory facility team, and has been 
observed since Chandra's 1999 launch. As a result, the Chandra Data Archive contains 2.16 million 
seconds of calibration observations acquired during 215 observing sessions (ObsIDs). Many of these 
observations were obtained to calibrate the off-axis performance of the detector. We use archive 
calibration data exclusively in this research, specifically those ObsIDs with high spatial resolution 
acquired on-axis between 2000 and 2009. With these combined data, we present new upper limits 
to the X-ray flux of a compact stellar object for the remnant. 



Observations 



We used Chandra archive ca libration data of E0 102 as observed with the Advanced CCD 
Imaging Spectrometer (ACIS; see iGarmire et al.l 2003, for a review of ACIS). Of the more than 
2 Ms of available Chandra calibration observations of E0102, we select those observations for which 
E0102 was imaged on the S3 chip in the ACIS-S conflguration. We also require that the calibration 
imaging was obtained without spectral gratings in the beam. This selection criteria ensured that the 
remnant was not imaged ac ross multiple chips, the optimal spatial r esolut ion of I'.'O was achieved 
with on-axis pointings (see IChandra Proposers' Observatory Guide! (120091 ) for more details), and 
also ensured a uniformity in the data set necessary for the spectral analysis presented in Section 4.2. 
The S3 chip is a back-illuminated CCD chip, and thus has higher sensitivity to lower energy photons 
than the front-illuminated 13 chip, making this chip the favorable one to use when characterizing 
any compact stellar object in the remnant with a soft spectral behavior (e.g. a Cas-A-like neutron 
star). The ObsIDs that met these criteria are listed in Table [TJ 



3. Data Reduction 

Data reduction and processing were completed with the Chandra software package, CIAO 
version 4.]|l|, and IDL. The ObsIDs were processed with the CIAO routine acis_process_events 
to remove pixel randomization. For each ObsID, a lightcurve was generated for a background 



^see |http : //cxc . hatrveLrd . edu/ciao/ 1 



region beyond the remnant, to check for possible flaring events. Temporal variations above the 
mean background sky, which could indicate solar flaring events during the exposure, were found 
to be minimal, for each of the ObsIDs selected. A charge transfer inefficiency (CTI) correction 
was applied, to minimize the effects of radiation contamination which the Chandra detectors have 
experienced in orbit since launch. Next, the sub-pixel event re-positioning algorithm ser_v2, an 
IDL procedure, was applied to each exposure. This routine is designed to improve C/i andra's spatia l 
resolution, by decomposing 2, 3, and 4 pixel photon splits at the chip interface (see iLi et al.l l2003l ) . 



Using the CIAO script merge_all, the data were re-projected onto one common frame (ObsID 
2843) and stacked to produce a single exposure with an effective exposure length of approximately 
324 ks. The increase in the signal-to-noise for this combined image is a factor of ~ 4 - 6 over the 
S/N of any individual ObsID. This stacked image is referred to as the merged image in the text. 
This, and all other, images which were used to complete the analysis we present here, at the native 
pixel scale of the ACIS detector 0'.'49 pixel'^ 



4. Analysis 
4.1. Image Analysis 

4.1.1. Locating a Compact Stellar Object 

Simulations and observations alike suggest that core-collapse supe rnova events are not likely 
to be spherically symmetric in their collapse (see IWang fc Wheeleiil2008l . and references therein). It 
may be possible, or even requisite for some varieties of core-collapse supernovae, that this asym- 
metric evolution r esults in a "pre-na tal kick" of the compact object revealed by the supernova event 



(JFryer et al.l 120061 : iFragos et al 



2OO9I ). We include this dynamic factor by defining an area in which 



a "kicked" stellar object would most likely be found. This "test source region" (hereafter, TSR) 



has a radius, r, estimated by r 
the remnant. 



a. 

te ' 



where Vk is the "kick velocity" and te is the expansion age of 



We use two published expansion ages to constrain te for E0102. From the analysis of the 
X-ray proper motions of shell ejecta over a period of ap proximately 20 year s, iHughes et al.l (l200d) 
determined the age of E0102 to be equal to ~ 1 kyr. lEriksen et al.l (I2OOII ) and iFinkelstein et al 



(J2006l ) determined from the analysis of medium and broadband Hubble Space Telescope imaging, 
that E0102 was significantly older, between ~ 2.0 and 2.1 kyr, respectively. These estimates to the 
age of the remnant agree within the error bars provided by each of the studies. To estimate the 
kic k velocity, V k , we a dopt the mean two-dimensional (2D) pulsar velocity of 240 km s^^, derived 
by iHobbs et al.l (|2005l ) from a proper-motion study of 73 pulsars with characteristic ages less than 
3Myr. With these assumptions, a "kicked" stellar remnant could travel a linear distance between 
0.1 and 0.5 pc. 

We assume the proper motion center (R.A. : 01^04'"02^05, decl. : -72°01'54'.'99 (J2000)) de- 



termined in iFinkelstein et al.1 (|2006l ) for the center of our TSR. The physical distance between this 
center and the main ejecta sheh of the remnant is approximately 6pc. Thus, a "kicked" stellar 
remnant with the aforementioned properties could not have reached the remnant shell in the ex- 
pansion lifetime. With the assumption of the IFinkelstein et al.1 (|2006l ) center, we define the TSR as 
a generous circular area with a radius ~ 1.2 pc, which we overlay on the map of the remnant in 
Figure [H 



4-1-2. Galactic Analogs to an E0102 Central Compact Stellar Object 

In the photometric analysis we perform here, we restrict our search for a stellar object in the 
rem nant to two well-studied potential analog s : the central compac t object (CCO) in Cassiopeia-A 



f see IChakrabarty et al. 



Crab pulsar (see 



2001 



Hester 



^ Fesen et al.ll2006l : iPavlov &: Lunall2009l . and references therein) and the 
20081 . and references therein). 



Gaetz et al.l (2000|) provide an upper limit to the X-ray luminosity of 9x 10^^ erg s ^ for a Crab- 



like pulsar in the remnant. This limit was determined from a 9.7 ks observation of th e remnant with 



the C handra ACIS instrument. Further, the spectroscopic study of this remnant by iFlanagan et al. 



(J2004l ) implies that if a Crab-like compact stellar object is present, then it must be faint at X-ray 
wavelengths. 

If the E0102 supernova event generated a Crab-like pulsar, then this object's spectral behavior 
would be well described by a power-law, with an approximately constant photon index (F) in the 

1.6 - 1.8, but the spectrum of the Crab remnant 



Chandra bandpass. For the Crab pul sar, F 



is softer and better fit with F ~ 2.0 (jWeisskopf et al.l 120041 : iSeward et al.l l2006l ) . At low energies 
(energy ~ 0.2 - 2.0 keV), the relative strength of a thermal emission component arising from the 
remnant plasma is likely to overwhelm the contribution from the Crab-like pulsar component. But, 
at higher energies, the relative strength of the Crab-like component would increase, possibly to the 
extent that the compact object could dominate the emission in the c enter of the remna nt, as was 



observed for the pulsar wind nebula associated with SNR G292+1.8 (JHughes et al. 



2001 



To identify such a trend for this remnant, we first processed the merged image to isolate high 
and low energy emission between 0.5 and 10.0 keV within energy bins whose widths were selected 
to include the following : known complexes of metal emission-lines of oxygen (O VII triplet : 0.46 
- 0.6 keV and O VIII Lya : 0.61 - 0.72 keV), neon (Ne IX triplet : 0.86 - 0.98 keV and Ne X Lya 
: 0.98 - l.lOkeV), magnesium (Mg XI triplet : 1.32 - 1.4keV and Mg XII Lya : 1.45 - 1.53 keV); 
those energies which provide upper and lower energy bounds for each of the previous line complexes 
(which could be used to identify continuum sources), i.e., a lower bound to oxygen ("Bin 1" : energy 
= 0.20-0.40 keV), a lower bound to neon and upper bound to oxygen ("Bin 2" : energy = 0.75-0.85) 
and a lower bound to magnesium and upper bound to neon ("Bin 3" : energy = 1.15-1.275 keV); 
and the hard (2.0-8.0 keV) continuum. We visually inspected the resulting X-ray spatial flux maps 
for each of these energy bins, but did not identify a compact object. A Crab-like pulsar should, in 



addition, be detected across a long wavelength baseline, but the aforementioned studies of E0102 in 
the radio, optical, and infrared wavelengths have not detected a neutron star. For these reasons, it 
is unlikely that a Crab-like neutron star is present in this remnant, and for the photometric analysis 
in this section we consider exclusively the presence of a Cas-A-like CCO. 



Chakrabarty et alJ (120011 ) combined ~ 6 ks of Chandra ACIS (S3 chip) observations of Cas-A, 
determined the location of the compact stellar object, and modeled the luminosity of a neutron star 
for the remnant. The energy of the maximum count rate for this object peaked at 1.6 ±0.2 keV 
and was well-fit with an ideal blackbody energy distribution (kT^g = 0.5 keV), assuming a column 
density A'^h = 1-1 x 10^'^ cm^^. The X-ray luminosity of the Cas-A CCO is equal to -^o.i-iokeV 
= 2.0lQ'g X lO'^^ erg s"^. Artificially moving this object to the distance of E0102 (59kpc), the 
unabsorbed flux in the same bandpass equals 4.7x10"^^ erg s~^ cm^, or equivalently a count rate 
equal to 7.15x10"^ counts s~^. In 324 ks, this count rate is equivalent to 221 counts. We cannot 
visually identify a compact object, with this count rate, within the TSR. 

An examination of the emission from the remnant within the TSR does reveal a strong gradient 
in the surface brightness, though. In Figure [H Region 1 contains an area within which the emission 
is of the lowest mean surface brightness within the TSR. In this 324 ks exposure, the mean count 
number within this class of region equals 117. The Q-spoke, the hot extended plasma feature 
extending from the proper motion center to the shell of the remnant, falls within the TSR and 
defines a "medium" class of surface brightness regions with average count rates ~ 5 times greater 
than the low surface brightness regions previously discussed. Regions 2 (mean count w 521) and 
3 (mean count ~ 459) are regions representative of this class of emission. Within a 2 x 2 square 
pixel region (center : R.A. - 0l'^04™02^015, deck- -72°01f55'.'82 (J2000)), the unique location of 
peak emission within the TSR defines a third class of surface brightness. The mean count value for 
this region. Region 4, equals ~ 1157 counts. This variation in surface brightness could imply that 
a point source exists within the remnant but that the remnant's hot plasma "masks" the emission 
that arises from the compact object. 

If we assume that a neutron star with a spectral behavior similar to the Cas-A CCO is in fact 
present in this remnant, but is concealed or masked by the X-ray bright plasma along the line of 
sight, then the observed count rates can be used to produce upper limits to the luminosity of such 
a neutron star. 



We assume a line of sight SMC HI Column Densit}!!, N^ = 5.82 x lO^i 



( 1999 ). With PIMMyfl, we then calculate the upper limit to the flux, and derive the luminosity, for a Cas-A 



neutron star with either a blackbody or power-law spectral energy distribution. For these calculations, we 
set the blackbody kTcH parameter and the power-law photon index F equal to the similar best-fit parameters 



^The Milky Way Column Density, estimated with the CXC software package GOLDEN an online tool accessible 
at |http : //cxc . harveird . edu/toolkit/colden . j sp[ is low relative to the SMC absorption, with TVh = 6.57 x 10^° 



'Available at |http : //cxc . heirvard . edu/toolkit/pimms . j sp] 



derived by IChakrabartv et al.l (|200l[ ). We present these upper limits in Table [2] In this table, the Scr rms 
sky noise (in counts) is never greater than '^ 3, which implies less than ^ 3% of the counts for each of the 
regions in Table [2] are detector, or sky, background not arising from the remnant. We do not attempt here 
to estimate the background "noise" , within each region, which arises from the remnant plasma directly (see 
Section 112]). 

A compact object is only one of many components contributing to total flux observed within a region 
of interest, and within the TSR (particularly within Region 4), the contribution from this source may be 
small relative to other components (background and sky "noise" , detector noise, detector response, absorbers 
along the line of sight, etc.). Nevertheless, for regions which are not coincident with the Q-spoke, the "zeroth- 
order" estimates to the upper limits to the flux from a compact object in these regions may be reasonable, 
as the contribution of a compact object in one of these regions cannot, by definition, exceed the observed 
flux. But, the TSR includes bright emission from the Q-spoke, and the surface brightness variations indicate 
that the contribution of the remnant plasma is variable with position. Thus, quantifying the contribution 
of a compact object to the total observed flux for each region requires a more rigorous estimation of the 
contribution of the background to the total observed flux. In the following sections, we present the results of 
two methods which provide a more robust determination of the relative strength of a compact object within 
the TSR. 



4- 1-3. Point-source Identification 

To date no study of E0102 has reported the detection of a compact object in physical association with 
the remnant. But, previous investigations of E0102 were limited in their sensitivity and spatial resolution. 
Combining the extensive archived observations of this remnant, we have produced a high signal-to- noise 
image of the remnant, at the superior X-ray spatial resolution of the Chandra ACIS instrument suite, which 
represents the best opportunity to identify a compact object within the TSR, particularly in the vicinity of 
Region 4. 

In the previous section, it was found that complex surface brightness variations (with radius, measured 
from the center of the remnant) prevent any conclusive visual identification of a compact object. In this 
section, we search for a point-source object in the image using the CIAO wavdetect script, which is the 
optimal source detection software package for Chandra X-ray images. This software iteratively identifies 
source objects by correlating image pixels with " Mexican Hat" wave lets, for a selection of user-defined 
wavelet scales(for a description of the software, see iFreemanet al.l 120021 ). We apply wavdetect for source 



detection in seven images of the remnant: Broad (energy = 0.5-7.0 keV), Medium (energy = 1.2-2.0 keV), 
Soft, (energy — 0.5-1.2keV), Hard (energy = 2.0-7.0 keVJj, and the three energy bands discussed in Section 
I4.1.2l that selected continuum emission between strong line complexes. For each of these images, wavdetect 
was used to identify sources on "Fine" (scales = "1 2") and "Coarse" (scales = "4 8 16") scale, with sigthresJ^ 
set equal to 4x10^^. We ran each wavdetect search twice, with and without an effective exposure map, 
and found the same number of detections within the TSR. 

WAVDETECT Searches in each of these images identified a source within Region 4, and the detection 

^These first four energy bands are directly comparable to Chandra Source Catalogue bands. 

^ sigthresh defines the threshold at which "cleaned" pixels are considered to be members of a source detection 



ellipse for each of these sources included part (or all) of this region. In Figure [2l we overplot the source 
detections which were identified within lO'.'O of the expansion center for the Broad, Hard, "Bin 2" and "Bin 
3" continuum images. In addition, we include in Figure [2] the "image file", an output from the wavdetect 
routine that is a representation of the data with background and noise subtracted. We are reluctant to 
conclude that these source detections are associated with a compact object, though, and believe that these 
results reinforce, quantitatively, the preliminary conclusion presented in Section [4.1.21 that a compact object 
may be present in or near to the brightest pixel region (Region 4) but masked, for the following reasons. 

First, the wavdetect results confirm that the Q-spoke is a real feature, contributes significantly to 
surface brightness variation across the center of the remnant, runs continuously between the center and 
outer edge (from the observer's perspective) of the remnant as illustrated in the "image file", and the 
brightest pixel region remains "masked" in the Q-spoke in the "image file" . 

WAVDETECT provides a variety of photometric measurements which can be used to constrain the ex- 
tended versus compact nature of detected objects (e.g., PSFRATIQj) and an estimate to the statistical 
significance of each detection (SRC_SIGNIFICANCE1J). When the source detections nearest to the brightest 
pixel region are considered individually, these statistical measurements could suggest possible point source 
detections, particularly in the narrow band images. Considered in the context of the full set of detections, 
though, the sources detected in the proximity of Region 4 are indistinguishable from other detections within 
the remnant, most of which are associated with hot gas background emission that could not be reduced from 
the image file by the wavdetect software. Furthermore, the lack of nearby point sources or double-peaked 
source detections make it difficult to interpret these results. In Table[3]we provide pertinent statistics for the 
set of detections from each of the wavdetect searches; here, we limit the set to those which have detection 
ellipse centers {a, 5) interior to the outer shell of the remnant. 

Lastly, each image's source detection ellipse most closely matching to, or including Region 4, in the 
"Fine"-scale wavdetect search is found to be offset from the center of Region 4 by less than I'.'O. But 
the size of the detection ellipse significantly increases at the "Coarse" scale. These latter two results likely 
indicate the difficulty which wavdetect has in adequately discriminating between sources and unresolved 
background components, or knots of hot plasma and "masked" compact objects. 

Though the WAVDETECT search results provide more point source candidates within the TSR, we believe 
that any photometric method of estimating the relative strength of a compact object "masked" by the Q- 
spoke is inadequate and now turn to a spectral analysis for characterizing the background emission. 



4.2. Spectral Analysis 

In the previous sections, we were limited in deriving firm upper limits on the detection and charac- 
terization of a compact object by the presence of a high surface brightness emission component that likely 



^WAVDETECT calculates this by dividing the estimated "radius" of the source cell to PSF^SIZE, an estimate of 
the point-spread function (PSF) calculated at the location of the detection. For further information, please see the 
CIAO Detect Reference Guide at |http://cxc. harvaLrd.edu/ciao/download/doc/detect_manual/wav_ref .htm l 

'^WAVDETECT provides an estimate to the significance of the detection with respect to the calculated background 
for the source detection region, but should not be taken at face value. This statistic is defined in the CIAO Detect 
Reference Guide and is calculated by dividing the net counts by the "Gehrels error" , ac = 1 + VBKGcovnts + 0.75. 



arises from the remnant plasma. In this section, we perform a spectral analysis designed to quantify the 
contribution of multiple background components present within an expanded set of regions within or near 
to the TSR. Only by fitting and reducing out the contribution of these components in the regions of interest 
will we be able to conclude firm upper limits for compact objects which may be masked by the Q-spoke. 

It is necessary to note that E0102 is used by the Chandra facility as a spectral calibration source, 
because it is a relatively constant and bright X-ray source, especially at energies between 0.55 and 0.7 keV 
corresponding to the emission-lines of O VII and O VIII. As such, it is an optimal target for u se in calibrating 



the re sponse of the ACIS CCDs and the variability of these CCDs response in time (see, e.g.. |Plucinskvetal 



(|2008[ )). Imaging of the remnant for calibration purposes does not preclude us from searching for multiple 



additional spectral components, though. 

Using the CIAO script SPECEXTRACT, we extracted spectra for each of the regions discussed in Section 
14.11 (see Figure [U and an additional six regions of varied size, shape, and location. We select these regions 
to better characterize the background emission within the remnant, not just within the Q-spoke and TSR. 
In Figure[3l we provide a map of the regions, interior to the remnant shell, for which spectra were extracted. 
Not displayed in Figure [3] is the background region (area ~ 27,000 pixeP) which was used to characterize 
the X-ray sky background not originating from the remnant itself. These regions were selected from areas 
of the chip where the individual ObsID exposures were, when stacked, spatially coincident at distance well 
beyond the remnant. For legibility. Figure |3] does not include Regions 2 and 3. 

In this section, we are concerned exclusively with the detection of a point source object by its spectral 
signature. In theory, the detection of any such object is possible only by identifying and reducing, from the 
spectral profile extracted for a particular region of interest, all emission which does not arise from the point 
source itself so that the flux from the compact object can be properly measured. We assumed that this 
background arose from one of three sources: the X-ray background sky arising from astrophysical sources 
independent of the remnant as well as instrumental background particular to the ACIS detector (i.e., "chip 
noise"), which we will refer to here as "Component A" background, and diffuse emission arising from E0102 
ejecta, the intensity of which varies with the position in the remnant, which we will refer to as "Component 
B" background. 

To calibrate the emission arising from the X-ray sky (e.g., diffuse galactic Milky Way and SMC emission), 
and the background arising from the ACIS detector electronics (e.g., the well- understood fluorescence lines 
of Si Ka, Al Ka, etc.; Baganoff 1999), we first analyzed a spectrum extracted from a uniform region (the 
aforementioned "background region") located approximately 1' from the outermost radial extent of the 
remnant. Using the XSPECO software package, we then modeled the broadband spectral behavior of this 
"Component A" emission using multiple power-law and Gaussian components to represent this diffuse X- 
ray sky emission and lines. With 583 degrees of freedom in the model fit to "Component A" background 
emission, the reduced minimized x^ statistic for the best-fit model is equal to 1.03. For this analysis, and 
throughout this section, we consider the spectral profile in two energy bins, "soft" X-rays (0.5 - 2.0 keV) 
and "hard" X-rays (2.0 - 8.0 keV ). 

"Component A" background emission is assumed to be uniformly present and of similar intensity 
within each of the regions for which we extracted spectra. This assumption does not hold for "Component 
B" emission which arises from the remnant itself. To fit this class of background emission, we must consider 



see http : //heasarc . gsf c . nasa . gov/docs/xanadu/xspec/ | 



each of the spectra individually. Because our goal is to identify a spectrally distinct compact object in 
this remnant, we aim to characterize the diffuse plasma spectrum rather than to physically understand it. 
Therefore, we fit an ad hoc model to the spectrum for each of the regions. Again, we used the xspec package 
to model this background emission, fitting the "Component B" X-ray emission (which was dominated, not 
unexpectedly, in the soft band by the strong oxygen and neon e mission-line comp lexes) with a system of 



Gaussian line profiles, and also with a variable APEC model fsee lSmithet al.ll200ll) . We acknowledge here 
that the spectral resolution is insufficient for distinguishing between spectral lines and continuum emission 
between those lines; thus the fits to these sources of background emission are model dependent. Therefore, 
in the following analysis, the upper limits to the continuum emission from any putative neutron star are 
implicitly model dependent. 

For each region, we combine both background models generated to produce a unified background model 
which we then fit to each of the observed spectra. In Figure SI we present the unified background model in its 
component form. In this figure, black data points with error bars represent the observed imaging spectrum 
extracted from, in this example, the brightest pixel region (Region 4). The "Component A" background 
model fit is overplotted in red: the red solid curve represents the sum of the individual model components 
(red dotted line) fit to the red observed data points (plotted with error bars) of the extracted spectrum from 
the background region. The unified background model, which includes the region-specific "Component B" 
background emission, is overplotted in black; in this figure the Component B emission model was generated 
using a variable APEC model and power-law components. For comparison, we provide in Figure [5] the 
background model fit, but with "Component B" emission-line complexes modeled with a system of Gaussian 
spectral lines. Each model reproduces the observations well, within the given errors, but in general the 
vAPEC models have smaller residuals. 

The best-fit background model was determined for each of the extracted spectra by iteratively mini- 
mizing the x^ statistic. The "Minimum x^" value for each best-fit unified background model is provided 
in Table ID Each model, for both the vAPEC and Gaussian fits, was defined with 577 degrees of freedom, 
thus the reduced x^ for each of the background fits was ^ 1.0 (e.g. Region 4, xt = 1.07). As "Compo- 
nent A" background emission was fixed in its contribution to the spectral energy distribution extracted for 
each region, this figure of merit essentially measures the goodness of the fit of "Component B" background 
emission, which arises predominantly from nebular line emission. 

Having sufficiently accounted for all sources of background emission, we can then derive an upper limit 
to the presence of a compact object within each of the regions. To derive these upper limits, we assume 
that a compact object is a neutron star and is well described by either a power-law or a blackbody model 
(hereafter, referred to as the "point-source model"). 

We define the point-source model as the product of a normalization constant and the specific functional 
form; for example, the power-law model point-source model function is Fe = /3 x E~ . Combining the 
best-fit unified background model, now fixed for each region by the previous analysis, with the point-source 
model, we define a two-component spectral model which we then re-fit to the observed spectrum for each 
region by only varying the magnitude of the normalization constant Q We increase the strength of the 
point-source component in this two-component model only until Ax^ = min.x^ + 6.63 is no longer satisfied. 
Note that 6.63 is the chi-square difference constant implying an observation is significant at the 99% (3 a) 
confidence level for a system modeled with 1 degree of freedom. Thus, with this analysis we have measured 



'We fix r (or kT, in the case of the Blackbody point source model) for each fit. 
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the upper limit to the observed flux from a compact stellar object, above which such an object would have 
been detected in this data. 

In Table |4] we present the results of this analysis, which we completed for three power-law models of 
a compact stellar object defined with F set equal to 1, 3, and 5, respect ively. This raii ge of indices spans 



the observable parameter space for the Cas-A CCO (power-law fit, see iPavlov fc Lunal '2009) and pulsars 
(i.e., the Crab pulsar). In Columns 2 - 6 of Table HI we present the flux, equivalent luminosity, and the 
minimum ^ statistic we determine from the best-fit unified background model to the spectrum extracted 
for each regions. In Columns 8-13 we present the upper limits to a neutron star for E0102. Specifically, 
in Columns 8 and 11, we present the flux difference between the two component best-fit model and the 
unified background fit for the region, for the soft and hard energy bands respectively. Positive values in 
both the soft and hard energy bands indicate that the normalization parameter that was found to satisfy the 
statistical criterion (defined previously) in one band, simultaneously satisfied the criterion in the other band. 
Or, similarly, the point-source model flux in both bands was sufficiently greater than the unified background 
model to be detected, minimally at 99% confidence. Negative or zero values in Columns 8 and 11 should be 
interpreted as non-detections at the 99% confidence level. These non-positive values can be attributed to 
a lack of discriminatory power in the models for differentiating between remnant hot gas emission and the 
point source component, within the observational errors, at the confidence level defined by the statistical 
criterion. The remnant is dominated by soft emission; thus these non-detections typically arise for the softest 
(r = 5) point source component model fits. We could produce upper limits for these models if we instead 
chose to enforce a stricter confidence level, e.g. 99.5%, which would require us to redefine the statistical 
criterion (in this example, Ax^ = min.x^ + 7.88). Lastly, in Columns 9 and 12 we present the flux upper 
limits, and in Columns 10 and 13 we present the equivalent luminosity upper limits for the stellar object. 

Because the brightest pixel region (Region 4 in Figure 2) has a distinct photometric signature (see 
Sections 14.1.21 and I4.1.3|) . we perform a similar analysis to that which was outlined above for this region, 
but also explore the possibility that the point source component within the two component model is best 
described by a blackbody model. We present this analysis for both the power-law (F = 1, 3, and 5) and 
Blackbody (kT^ff — 0.25, 0.5, and 1.0) model point source components in Table [5l 



5. Discussion 

We interpret the non-detection of a compact stellar object in the remnant E0102 within the context 
of three physical scenarios. In the first of these scenarios, the non-detection of a compact stellar object for 
E0102 implies that such an object simply does not exist. There are two plausible scenarios to consider when 
addressing this possibility. 

First, a stellar compact object may have been destroyed during the supernova event. If E0102 repre- 
sents the e nd stage of a Type la supernova event, then the progenitor white dwarf (s) would have likely been 



destroyed (IWooslevetal.lll986f) . But, it is difficul t to reconcile this formation scenario with the evidence pre- 



sented bv lBlairet al.l ( 2000 ) and lChevalieiJ ( 20051) which would support a Type Ib/c or Ilb/L class supernova 



event for the E0102 progenitor. 

The theory of Type Ib/c and Type II supernova predicts that compact objects are generated during 
the supernova event, and compact stellar objects have been observed in physical association with ~ 25% of 
Galactic remnants (see Green's Catalogue] ?). But, in addition to these neutron-degenerate stellar objects. 
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core- collapse supernovae of sufficiently massive (> 20Mq) progenitors may generate black holes (see ISmartt 
20091 . for a review). At least two arguments exist in support of the black hole sce nario, and though the; 
are circumstantial they should be addressed in future work on the remnant. First, iFinkelstein et al.l ( 200' 
suggest that E0102's progenitor star may have been a massive (>50 Mq) Wolf-Rayet star, based on the 
consideration of th e local environment o f the remnant. Second, to explain a similar remnant morphology 
in the SNR W49B 



model; 



Keohane et al.l (|2007|) suggest a possible jet-driven explosion mechanism (a "collapsar" 



McFadvenet al.ll2000l ). which would generate a black hole. 



A second interpretation of the non-detection of a compact stellar object in the remnant is that a neutron 
star in this remnant is in fact present, but is both faint in its emission and obscured by the Q-spoke. This 



possibility is motivated by the imaging and analysis presented in Sections 14.1.21 14.1.31 and 14.21 

A final interpretation of the non-detection of a compact stellar object is systemic, such that the estimate 
to Njj is too low. Thus, the flux from a compact object in this remnant is absorbed by the interstellar medium 
(ISM) along the line of sight. At a distance of 59kpc, the Cas-A CCO would have a detectable signal; in this 
324ks ACIS image, this signal would be equivalent to ~ 240 counts between the energies of 0.1 and lOkeV. 
Such a source would stand well above the noise in this observation. But, the average number count within 
that portion of the TSR unobscured by the Q-spoke is equal to ~ 140 in the 324 ks exposure. By measuring 
the relative strength of the background in an image filtered for energies greater than ~ 10 keV (above which 
the effective area of Chandra + ACIS is nominally zero with respect to the known background for the S3 
chip), an estimate to the cosmic ray background count rate can be made. With this rate appropriately scaled 
to the exposure length and area of this region, it can be shown that only a small fraction of these counts (<< 
1 %) can be attributed to cosmic rays. Thus, to reduce the Cas-A-like neutron star fiux to this observed 
level would require a column density along the line of sight ^12 times higher than we assumed. This is not 
reasonable given the errors on Njj , and we conclude that gas along the line of sight is not likely to obscure 
the flux from a compact object in this remnant. 

The non-detection of a compact object for this remnant, though, does not preclude its detection at 
shorter or longer wavelengths than the Chandra operational energy bands. In particular, this research 
suggests that imaging of the remnant at th ose wavelengths at w hich emission from cjccta material is relatively 
low may reveal such an object. Recently, ICamiloet al.l (|2009l ) presented the detection of a radio pulsar in 
association with Galactic SNR G315. 9-0.0. This remnant displays a similar emission "spoke," extending 
perpendicularly outward from the remnant, with the radio pulsar detected at the outermost tip of Q-spoke. 
The remnant E0102 and G315. 9-0.0 should not be directly compa red because of morphological differences 
and different estimated ages, but the results of lCamiloet al.l (J2009|) and this research suggests that a more 
extensive review at radio wavelengths could provide insight into the nature of a compact object for this 
remnant. Furthermore, the recent result of 16 galactic pulsars discovered with the Fermi Gamma Ray 
Telescope (| Abdo et al.l |2009|) suggests that imaging of the remnant at higher energies than the Chandra 
bands may be fruitful. In addition, a high- resolution timing study of the center of the TSR could detect 
a pulsar. At present, Chandra GO program 11500201 (PI : Petre) is approved for 20 ks of imaging with 
HRC-I in Spring 2010. These data could be used in the assessment of temporal variability of an, as of yet, 
undetected central pulsar. 
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6. Conclusion 

We have combined 324 ks of Chandra X-ray observations of the SNR 1E0102.2— 7219 in the SMC to 
produce the highest signal-to-noise image of this young, oxygen-rich remnant. In this research, we attempt 
to identify a compact stellar object for this remnant within a generous search area motivated by observations 
of this remnant. Galactic SNRs, and compact stellar objects. We cannot confirm the presence of a compact 
object for this remnant, but are able to provide the strongest upper limits on the flux, assuming that the 
spectral behavior of a faint neutron star for this remnant is best described by a blackbody or power-law 
spectrum energy distribution. 
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Fig. 1. — Remnant 1E0102.2-7219, displayed at the native pixel scale of the ACIS instrument (0'.'49 
pixel"^). Emphasized here is the region in which a compact stellar object would likely be found 
(see Section [4.1.ip . The upper limits to the flux from a neutron star in each of these regions are 
discussed in Section 14.1.21 
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Table 1. E0102 Calibration Observations with S3 Chip Events Providing at Least l"Resolution 



Year 


Exposure Time 


(ks) 


ObsIDs 


2000 


35.00 




1308, 1423, 1803 


2001 


32.48 




141, 1311, 2843, 2844 


2002 


15.33 




2850, 2851 


2003 


51.70 




3519, 3544, 5123, 5251, 5252 


2004 


55.12 




5130, 5131, 6074, 6075 


2005 


26.75 




6042, 6043 


2006 


53.30 




6758, 6759, 6765, 6766 


2007 


20.98 




8365 


2008 


19.20 




9694 


2009 


14.57 




10655, 10656 



Table 2. Upper Limits to the Flux of a Cas-A-like Neutron Star in Remnant 1E0102. 2-7219 



Region^ 


Count Rate^ 


Flux'' 


Luminosity "* 


-Lx 5 


Flux 


Luminosity 


Lx 




3.61 (-4) 


Assume T =2.35 




Assume F =3.13 




1 


1.12 (-14) 


4.65 (33) 


0.62 


3.73 (-14) 


1.55 (34) 


0.35 


2 


1.61 (-3) 


5.01 (-14) 


2.09 (34) 


2.81 


1.66 (-13) 


6.89 (34) 


1.60 


3 


1.41 (-3) 


4.39 (-14) 


1.82 (34) 


2.46 


1.46 (-13) 


6.06 (34) 


1.41 


4 


3.57 (-3) 


1.11 (-13) 


4.61 (34) 


6.22 


3.69 (-13) 


1.53 (35) 


3.56 




3.61 (-4) 


Assume kT=0.53 




Assume kT=0.49 




1 


3.49 (-15) 


1.45 (32) 


0.85 


3.46 (-15) 


1.44 (33) 


0.72 


2 


1.61 (-3) 


1.56 (-14) 


6.47 (33) 


3.81 


1.55 (-14) 


6.43 (33) 


3.21 


3 


1.41 (-3) 


1.37 (-14) 


5.69 (33) 


3.34 


1.35 (-14) 


5.60 (33) 


2.80 


4 


3.57 (-3) 


3.46 (-14) 


1.44 (34) 


8.45 


3.43 (-14) 


1.42 (34) 


7.12 



Note. — (-k) is equivalent to x 10 ^ 

^For the location of each region, see Figure [l] '^ in counts s~ 



Unabsorbed (0.1-10 keV) Flux 



[erg s ^ cm ^] 



derived from Observed count rate, and calculated by PIMMS assuming Cycle 10 
systematics for ACIS-S with no gratings; * Unabsorbed Luminosity [erg s~^]; ® Lx ■ Estimate of 
the Unabsorbed Lumin osity from previous column, Lc : Unabsorbed Luminosity determined by 
Chakrabartvetal.1 (120011 1. 
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-l.Spc 



Fig. 2. — Displayed (from top to bottom left hand column) are the broad and hard broadband 
images and "Bin 2" and "Bin 3" narrowband images of E0102 at the ACIS native pixel scale 
(0'.'49 pixel"^) in gray scale. The online edition of this figure is provided in color. Overplotted on 
each image are the sources detected within IC/.' of Region 4 (see Figure 3) with CIAO package 
WAVDETECT for the "Fine" -scale (in green) and "Coarse" -scale (in blue) wavdetect searches. A 
red region is also plotted to indicate the center of Region 4. Displayed in the right column is the 
image file produced by wavdetect which represents the remnant with the estimate of background 
and sky noise subtracted. Definitions of the energy bands on which these images were binned the 
WAVDETECT searches performed and the results and discussion derived from these source detections 
can be found in Sections 4.1.2 and 4.1.3. 
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~1.5pc 



Fig. 3. — Overlaid on this image of the remnant are those regions for which a spectral analysis was 
performed; we refer the reader to Section 14.21 for complete details and results of this analysis. The 
remnant is displayed at the native pixel scale of the ACIS instrument ((y.'49 pix~^). 



Table 3. Source Detections with wavdetect in E0102 



Search Scale 



Significance 



PSF Ratio" Min. PSF Ratio' Signiflcance" Max. Signiflcance" 



Nearest Detection 



Remnant Detections 



Coarse 
Fine 



5.72 
1.98 



176.4 
98.5 



11 
93 



6.9 ± 2.1 
2.2 ± 0.8 



3.41 
1.6 



544 ± 259 
31 ± 32 



Sec S ect ions 14.1.21 and f^. 1.3 1 for energy band definitions. 
Fine" = WAVDETECT "scale" parameter — "1 2"; "Coarse" = WAVDETECT parameter "scale" — "4 8 12". 

WAVDETECT parameter PSFRATIO which can be used identifying extended sources; sec Scction l4.1.3l 
*WAVDETECT parameter SRCSIGNIFICANCE; see Section ITTTsl 

Number of detections within ^ 25. of the E0102 expansion center. 

Mean, and la variance from mean, of WA\'DETECT parameter PSFRATIO for all detections (N ^Nre-m) in energy band. 
^Minimum PSFRATIO value for source detection interior to E0102. 
^Maximum SRCSIGNIFICANCE value for source detection interior to E0102. 



922 
132 



Soft 


Coarse 


5.73 


159.9 


10 


7.1 ± 1.8 


3.46 


556 ± 212 


887 


Fine 


1.95 


48.2 


91 


2.1 ± 0.8 


1.4 


31 ± 32 


139 


Hard 


Coarse 


5.56 


16.4 


15 


6.0 ± 2.3 


5.10 


38 ± 19 


76 


Fine 


1.95 


8.1 


31 


2.2 ± 0.7 


1.75 


6 ± 2 


10 


Medium 


Coarse 


6.85 


130.4 


13 


6.9 ± 2.0 


3.42 


172 ± 74 


302 


Fine 


3.21 


36.7 


76 


2.4 ± 0.8 


1.8 


15 ± 10 


45 


"Bin 1" 


Coarse 


7.17 


11.5 


11 


7.1 ± 2.01 


3.09 


51 ± 17 


79 


Fine 


1.96 


6.4 


26 


1.8 ± 0.9 


1.8 


6 ± 3 


14 


"Bin 2" 


Coarse 


3.11 


47.7 


10 


7.2 ± 2.3 


3.11 


187 ± 71 


308 


Fine 


3.13 


29.6 


67 


2.2 ± 0.8 


1.8 


15 ± 11 


46 


"Bin 3" 


Coarse 


6.83 


44.4 


11 


6.9 ± 1.9 


4.74 


107 ± 47 


172 


Fine 


3.20 


23.3 


63 


2.1 ± 0.9 


1.4 


10 ± 6 


30 



CO 

I 



Detection nearest, in proximity, to Region 4. 



Table 4. Upper Flux Limits from Spectral Analysis to an E0102 point source: Power Law Model 



Region 



Mill, xl * 



r, Photon Indcx^^ 











3ackgrou 


iid Fit 














Point 


Source 


Fit 




























1 


9.600 ( 


-16) 


2.263 


-14) 


9.391 


33) 


3.800 (-15) 


2.471 


-13) 


1.025 


(35) 


pi 


2.167 


-14) 


8.992 (33) 


2.433 


(-13) 


1.010 (35) 


1.07 


3 


9.000 ( 


-16) 


2.257 


-14) 


9.366 


33) 


2.000 (-16) 


2.435 


-13) 


1.010 


(35) 


















5 


3.600 ( 


-16) 


2.203 


-14) 


9.142 


33) 


0.000 


2.433 


-13) 


1.010 


(35) 


















1 


7.900 ( 


-16) 


2.281 


-14) 


9.465 


33) 


3.100 (-15) 


2.438 


-13) 


1.011 


(35) 


pS 


2.202 


-14) 


9.138 (33) 


2.407 


(-13) 


9.988 (34) 


1.05 


3 


7.300 ( 


-16) 


2.275 


-14) 


9.441 


33) 


1.000 (-16) 


2.408 


-13) 


9.992 


(34) 


















S 


-1.500 


(-16) 


2.187 


-14) 


9.075 


33) 


8.900 (-15) 


2.496 


-13) 


1.036 


(35) 


















1 


8.200 ( 


-16) 


2.267 


-14) 


9.407 


33) 


3.300 (-15) 


2.476 


-13) 


1.027 


(35) 


rl 


2.185 


-14) 


9.067 (33) 


2.443 


(-13) 


1.014 (35) 


1.08 


3 


9.000 ( 


-16) 


2.275 


-14) 


9.441 


33) 


2.000 (-16) 


2.445 


-13) 


1.015 


(35) 


















5 


3.600 ( 


-16) 


2.221 


-14) 


9.216 


33) 


0.000 


2.443 


-13) 


1.014 


(35) 


















1 


8.000 ( 


-16) 


2.302 


-14) 


9.552 


33) 


3.100 (-15) 


2.488 


-13) 


1.032 


(35) 


r2 


2.222 


-14) 


9.221 (33) 


2.457 


(-13) 


1.020 (35) 


1.05 


3 


7.900 ( 


-16) 


2.301 


-14) 


9.548 


33) 


1.000 (-16) 


2.458 


-13) 


1.019 


(35) 


















5 


3.600 ( 


-16) 


2.258 


-14) 


9.369 


33) 


-1.000 (-16) 


2.456 


-13) 


1.019 


(35) 


















1 


8.000 ( 


-16) 


2.302 


-14) 


9.553 


33) 


3.200 (-15) 


2.389 


-13) 


9.913 


(34) 


rS 


2.222 


-14) 


9.221 (33) 


2.357 


(-13) 


9.781 (34) 


1.06 


3 


9.000 ( 


-16) 


2.312 


-14) 


9.594 


33) 


2.000 (-16) 


2.359 


-13) 


9.789 


(34) 


















5 


4.200 ( 


-16) 


2.264 


-14) 


9.395 


33) 


0.00 


2.357 


-13) 


9.780 


(34) 


















1 


8.000 ( 


-16) 


2.300 


-14) 


9.544 


33) 


3.200 (-15) 


2.467 


-13) 


1.023 


(35) 


r5 


2.220 


-14) 


9.212 (33) 


2.435 


(-13) 


1.010 (35) 


1.05 


3 


8.900 ( 


-16) 


2.309 


-14) 


9.581 


33) 


2.000 (-16) 


2.437 


-13) 


1.011 


(35) 


















5 


3.000 ( 


-16) 


2.250 


-14) 


9.336 


33) 


0.000 


2.435 


-13) 


1.010 


(35) 


















1 


8.100 ( 


-16) 


2.294 


-14) 


9.519 


33) 


3.200 (-15) 


2.476 


-13) 


1.027 


(35) 


re 


2.213 


-14) 


9.183 (33) 


2.444 


(-13) 


1.014 (35) 


1.05 


3 


8.300 ( 


-16) 


2.296 


-14) 


9.527 


33) 


2.000 (-16) 


2.446 


-13) 


1.015 


(35) 


















5 


3.600 ( 


-16) 


2.249 


-14) 


9.332 


33) 


0.000 


2.444 


-13) 


1.014 


(35) 


















1 


8.100 ( 


-16) 


2.281 


-14) 


9.465 


33) 


-7.000 (-15) 


2.269 


-13) 


9.416 


(34) 


rl 


2.200 


-14) 


9.129 (33) 


2.339 


(-13) 


9.706 (34) 


1.05 


3 


8.400 ( 


-16) 


2.284 


-14) 


9.477 


33) 


-1.000 (-14) 


2.239 


-13) 


9.291 


(34) 


















5 


3.900 ( 


-16) 


2.239 


-14) 


9.291 


33) 


-1.030 (-14) 


2.236 


-13) 


9.278 


(34) 



o 



Note. — ( - k) is equivalent to X 10 

T:Modeled background flux "s" designates the soft X-ray band { 0.5 — 2.0keV ) and "h" designates the hard X-ray band ( 2.0 — S.OkeV ) modeled flux (erg s~ cm~"") — ^T .;gquiYalent 
background luminosity, assuming distance to E0102 equal to 59 kpc; "s" designates the soft X-ray band { 0.5 — 2.0 keV ) and "h" designates the hard X-ray band ( 2.0 — S.OkeV ) modeled flux 
(erg s~ cm~ ) — ^ :Minimum Xj^i i-'^577; please see Section |4.2| — "" : Assumed Power Law Index, P — : the diff^erence between the modeled unified background flux and the flux of the model 
point source, in energy band x — '^ the point source upper flux limit, in energy band x — : the equivalent point source upper luminosity, in energy band x 



Table 5. Upper Limit Flux for the Brightest Pixel Region' in E0102 



Min. x^** 



Fixed Parameter^ 



Background Fit 



Point Source Fit 



Blackbody Fit 2.171 (-14) 9.009 (33) 2.309 (-13) 9.581 (34) 0.78 



0.25 0.163 (-14) 2.334 (-14) 9.685 (33) 0.001 (-13) 

0.5 0.113 (-14) 2.284 (-14) 9.478 (33) 0.008 (-13) 

1.0 0.061 (-14) 2.232 (-14) 9.262 (33) 0.025 (-13) 



2.310 (-13) 9.585 (34) 
2.317 (-13) 9.615 (34) 
2.334 (-13) 9.685 (34) 



Power-Law Fit 2.199 (-14) 9.125 (33) 2.442 (-13) 1.010 (35) 1.07 



7.900 (-16) 2.278 (-14) 9.453 (33) 3.100 (-15) 
7.600 (-15) 2.959 (-14) 1.228 (34) 1.00 (-16) 

3.300 (-16) 2.232 (-14) 9.262 (33) 0.00 



2.473 (-13) 1.026 (35^ 
2.443 (-13) 1.014 (35) 
2.442 (-13) 1.013 (35) 



Note. — ( — k ) is equivalent to X 10 



t 



Region 4, sec Figure^ — ^T .]y[Q(j(-,iQ(^ Background Flux "s" designates the soft X-ray band ( 0.5 — 2.0 keV ) and "h" designates the hard X-ray band ( 2.0 — 8.0 keV ) modeled flux {erg t 



')- 



*:Equivalent Background Luminosity, assuming distance to E0102 equal to 59 kpc; "s" designates the soft X-ray band ( 0.5 — 2.0 keV ) and "h" designates the hard X-ray band ( 2.0 — 8.0 keV ) modeled 
flux (erg s~ cm~-) — ^^:Minimum x^i ^' — 577; please sec Scction |4.2| — "^lAssumcd model parameter, F or kT fx (keV) — ■.A.Fx is equal to difference between the modeled unified background flux 
and the flux of the model point source in energy band x — '^: the point source upper flux limit, in energy band x — : the point source upper luminosity liiTiit, in energy band x 
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Fig. 4. — "Soft" imaging spectrum extracted from the brightest pixel region in the remnant (Region 
4, see Figure [3j) is plotted with observational error bars in black. Overplotted (in red and black) 
are the best-fit models to the background emission arising from all sources. Here, the nebular line 
emission complexes ("Component B" emission) have been modeled with a variable APEC thermal 
plasma model in xspec. Please see Section [4.21 for details on the sources of background emission 
and the results of the spectral analysis performed for each of regions identified in Figures [1] and [3j 



23 



> 

0) 



U) 



c 

Z5 

o 
o 

"O 

CD 
N 

O 

E 

i^ 

o 

c 



D 

zs 

CD 



0.015 



0.01 



5x10" 



2x10 





-3 





-2x10" 



f 1 


Mm! 






III 




1 1 

j . '" 1 ' 1 'MM 1 ' 


::;:-V-:-" r-^ -^f-^'"-- -^^ ■ ■- 1- -:---::-'— '-^^"-'-y:^.---v\m m 


wW™[^i'hi'v,;,;iJ,',,,ii'irf:,i,i,iliM^ 












IlllJ 


J!,li, lllill|ili<l|iiii!i<iii<ii|iiiiiiiiliiliili 1 III 


1 1 1 


1 -- 1 1 


tfl'flllll 




|||.^i.||Jii|Tv™iii;,';iJi;';,vi,|i,i((|i|;i,!4f^ 

1 II" 



0.5 



2 



Energy (keV) 



Fig. 5. — Same as Figured! the "soft" imaging spectrum extracted from the brightest pixel 
the remnant (Region 4, see Figure [3]) is provided here, with model fits to the background 
overplotted. Here, the hot gas background component arising from the remnant has been 
with Gaussian spectral lines in xspec. Please see Section [4.21 for more details. 
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